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ABSTRACT 

Next generation air and space vehicle designs are being 
driven by increased performance requirements, demanding a high 
level of design integration between traditionally separate design 
disciplines. Interdisciplinary analysis capabilities have been 
developed for aeroservoelastic aircraft and large flexible 
spacecraft , but the requisite integrated design methods are only 
beginning to be developed. One integrated design method which 
has received attention is based on hierarchal problem 
decompositions, optimization, and design sensitivity analyses. 
This paper highlights a design sensitivity analysis method for 
Linear Quadratic Gaussian (LQG) optimal control laws, enabling 
the use of LQG techniques in the hierarchal design methodology. 
The LQG sensitivity analysis method calculates the change in the 
optimal control law and resulting controlled system responses 
due to changes in fixed design integration parameters using 
analytical sensitivity equations. Numerical results of a LQG 
design sensitivity analysis for a realistic aeroservoelastic aircraft 
example are presented. In this example, the sensitivity of the 
optimal control law and aircraft response for various parameters 
such as wing bending natural frequency is determined. The 
sensitivity results computed from the analytical expressions are 
used to estimate changes in response resulting from changes in 
the parameters. Comparisons of the estimates with exact 
calculated responses show they are reasonably accurate for ± 
1 5% changes in the parameters. It is also shown that evaluation 
of the analytical expressions is computationally faster than 
equivalent finite difference calculations. 

INTRODUCTION 

The design of new generation air and space vehicles is 
increasingly becoming subject to extensive requirements for 
design integration, that is, the close coordination of the design of 
the various pans of the vehicle. For example, many modem 
fighter aircraft require integration of the flight controls and 
engines so that sufficient power is available at all flight 
conditions. Similarly, the aircraft flight control and structural 
designs must be integrated to avoid potential aeroservoelastic 
instabilities. To meet the challenge of integrated aircraft design 
requirements, methods which tie together existing aerodynamic, 
structure, control, and propulsion design methods are needed. 

One such integrated design methodology currently under 
development at the NASA Langley Research Center is based on 
hierarchical problem decompositions, multilevel optimization 
methods, and design sensitivity analyses. 1 This methodology 
depends on the decomposition of the integrated vehicle design 
problem into vehicle requirements and separate aerodynamic, 
structure, control, and/or propulsion subsystem requirements. 
The subsystem designs are obtained independently subject to a 
set of fixed design integration parameters, using existing design 
methods and tools. An iterative optimization method is used to 
satisfy the integrated vehicle design requirements through 
modification of the design integration parameters and repeated 
subsystem designs. Subsystem design sensitivity data relative to 
the design integration parameters are used as the gradient 
information for the optimization procedure. 

An application of the hierarchal integrated design 
methodology is to the aeroservoelastic design of aircraft control 
laws and structure, including the effects of unsteady aerodynamic 
forces due to structural and control surface motions. This 
application requires the use of dynamic response requirements as 
the integrated design objective and a control law design method 


that uses the feedback signals actually available from the aircraft 
sensors. Both of these requirements necessitated the 
development and validation of appropriate design sensitivity 
expressions. Linear Quadratic Gaussian (LQG) optimal control 
law methods were selected for the control law design. Dynamic 
response criteria considered include time responses to control 
surface motions and discrete aerodynamic gusts, stochastic 
responses to random gust environments, closed-loop system 
eigenvalues, and open- and closed- loop frequency responses. 

The sensitivity developments have recently been 
completed. 2 A summary of these results and the application and 
validation of the sensitivity expressions to an aeroservoelastic 
aircraft example are described in this paper. Sensitivity results 
have been computed and are shown for design integration 
parameters related to aircraft wing bending stiffness, feedback 
accelerometer location, and control law design specifications. 
TTiese parameters are typical of those which would be used to 
obtain an integrated structure/control law design of an 
aeroservoelastic aircraft by the hierarchal design method. The 
sensitivity results were also used to validate the analytical 
sensitivity expressions. This was accomplished by comparing 
the sensitivity result with changes in responses due to control 
laws designed for different fixed values of the design integration 
parameters. Finally, the relative computational effort of 
computing the sensitivity information using analytical 
expressions versus numerical finite difference methods was 
investigated. 

INTEGRATED DESIGN METHOD 

A general integrated structure/control law design 
formulation based on hierarchal decompositions and multilevel 
optimization is shown in Figure L In Figure 1, the structural 
design and the control law design are independent lower level 
design problems. These lower level designs are coordinated 
using a set of design integration parameters. The upper level 
design optimization problem reflects the desired objectives of the 
integrated aircraft smicture/comrol law design. As a hypothetical 
example, the upper level objective might be to reduce peak 
transient responses of the aircraft due to a gust encounter and to 
reduce the weight of the structure. The actual peak transient 
responses of the aircraft would come from analysis of the control 
law design at the lower level, while the actual structural weight 
would come from the lower level structural optimization. These 
might then be combined as a weighted sum of square errors 
between the actual and desired values to form a single upper level 
optimization performance index. The upper level design 
variables, which are the design integration parameters, would 
then be selected to optimize the integrated design. 

The values of the design integration parameters at any 
time are treated as fixed for the lower level designs. The 
sensitivities of the lower level designs to these fixed parameters 
are computed and used in turn to compute the gradient of the 
related part of the upper level performance index. That is, these 
sensitivities are the gradients necessary to perform the top level 
optimization. In the present hypothetical example, one of the 
design integration parameters may be a local structural stiffness 
requirement, which appears as an equality constraint in the lower 
level structural design. The sensitivity of the optimized structural 
weight to this parameter is computed at the lower level and 
returned for use in computing the part of the gradient of the upper 
level performance index that is related to structural weight. 
Another of the design integration parameters might be a 
maximum allowable control surface deflection limit. The 
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sensitivity of the optimal control law design with respect to this 
parameter would then be used to compute the sensitivity of the 
peak transient gust response of the controlled aircraft, as required 
to perform the upper level design integration optimization. 

In many cases, existing nonlinear programming- based 
structural optimization and design sensitivity analysis methods 
can be used for the lower level structural design. These methods 
may themselves be hierarchal, multilevel optimization 
algorithms. 3 ^ 


In the rest of this paper, the use of Linear Quadratic 
Gaussian optimal control law design methods in hierarchal 
integrated aircraft structure/control law design is examined in 
detail. Expressions for the sensitivity of controlled system time, 
frequency, and stochastic responses in terms of state-space 
coefficient sensitivity matrices are discussed below. The 
sensitivity of optimized LQG control law to fixed parameters 
must be known to compute the needed state-space coefficient 
sensitivity matrices. Analytical expressions for the sensitivities 
of the LQG gain matrices to fixed problem parameters are 
discussed next, followed by the controlled system response 
sensitivity expressions, 

LQG CONTROL LAW SENSITIVITY 


Analytical expressions for the sensitivity of the optimal 
LQG problem solution above to the parameters pi can be obtained 
by differentiation of the necessary conditions of optimality. 
Since the solution of the LQG optimal control problem is the 
interconnection of the optimal Linear Quadratic Regulator (LQR) 
and Kalman Filter (KF), the necessary conditions for the LQG 
problem are the necessary conditions for the LQR and KF 
problems. Detailed derivations for the sensitivity expressions are 
presented in reference 2. The results arc summarized here. 


Regulator Sensitivity - The sensitivity of the optimal 
LQR gain matrix G s with respect to the 1 th parameter pi is 

= . R '1R r ' B y S + R '^i S + R ' (6) 

op dp dp 1 dp 

where the subscript i is dropped for convenience throughout the 
remainder of the paper except where necessary to avoid 
confusion. The unknown sensitivity of the steady-state LQR 
Ricatti solution S is obtained from 
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The Linear Quadratic Gaussian (LQG) optimal control 
law problem 5-7 is to find the control u(t) for the system 
x s 38 A s Xj + B s u + D s w s (la) 

y = C,x s (lb) 

z = M s x s + v (lc) 

such that the cost function 

T 

J = E j[y T Qy + u T Ru) dt (2) 

o 

is minimized, where E denotes expected value and x is the final 
time. In equations ( 1 ), x s is the system state vector of dimension 
(n, x 1), y is vector of system responses of dimension (n y x 1), 
and z is a vector of measured system outputs of dimension (n 2 x 
l). The vectors wj (n w x 1) and v (n v x 1, n v = n z ) are zero 
mean, Gaussian distributed, white noise disturbances with 
intensity matrices W s and V respectively, and the matrices A s , 
B s , C s , D s , and M s are real valued coefficient matrices of 
appropriate dimensions. 

It is assumed that the matrix pair (A s , B*) is stabilizable, 
the pair (M s , A s ) is detectable, and the pair (Q 0 . A s ) is detectable, 
where Qo is defined by Q 0 T Qo = C S T QQ and the matrices Q and 
R are positive semi-definite and positive definite respectively. 
The solution of the LQG optimal control law problem is then the 
interconnection of the optimal Linear Quadratic Regulator and the 
optimal linear state estimator or Kalman Filter. 5-7 In the Kalman 
Filter, the measured outputs z arc used to create estimates for the 
actual system states x s . Thus, the optimal control law is 


u = - G s x f (3a) 

x f = A s x f + B s u + Fj(z - M s x f ) (3b) 

where the gain matrices Ff and G s are given by 
F ( = TM*V' (4a) 

G, = R 'bJs (4b) 

and the matrices S and T are the positive definite solutions of the 
steady-state nonlinear matrix Riccati equations 

0 = aIs + SA, - SB.r 'bIs + cloc, (5a) 

0 = A S T + Ta[ - Tm[v‘m,T + D.W.dJ" (5b) 


It is assumed that the state-space model coefficient 
matrices A s . B s , C s , D s , and M s , the noise intensity matrices W s , 
and V, and the cost function weighting mainces Q and R arc 
time- invariant, continuous, differentiable functions of a number 
of parameters pj, i = l,...,nn, whose nominal values arc fixed 
during the solution of the LQG optimal control law problem. It is 
further assumed that the functional dependence of the above 
matrices on the parameters is known so that the partial derivative 
of each matrix with respect to each parameter i* ^ ^ known. 


Equation (7) is a linear Lyapunov equation which has a 
unique solution [reference 5, pg. 103, Lemma 1.5] since the 
coefficient matrix A s - B S G S is asymptotically stable by the 
properties of the LQR solution [reference 5, pg. 237, Theorem 
3.7]. This equation must be solved for each different design 
parameter pj, however the coefficient matrix (A s - B S G S ) is the 
same for every parameter. This can be used to advantage in 
developing a numerical algorithm for solving equation (7) for 
large numbers of parameters. 


Kalman Filter Sensitivity - The sensitivity of the 
Kalman Filter gain matrix Ff with respect to p is 

+ T^-V' 1 - TmJv 'I^V 1 (8) 

dp dp s t)p dp 

The unknown sensitivity of steady-state KF Ricatti solution T is 
obtained as the solution of 



+ - mIv l gv 1 M s + } (9) 

which is also a linear Lyapunov equation with a unique solution 
since the coefficient matrix A s - FfM s is asymptotically stable by 
the properties of the Kalman Filter solution. The coefficient 
matrix is the same for every parameter pi for this Lyapunov 
equation as well. 


Optimal Controlled System Sensitivity - Defining 
a state estimate error vector e, an augmented state vector x, and 
an augmented noise vector w as 
e = x s - x f 


(10a) 
(10b) 
(10c) 

the open- loop combined system and state estimator can be written 
in state-space form as 

i = A 0 x + Bu + Dw (1 la) 

y = Cx Gib) 


J*.l 

-PI 


u = Gx 

where the matrices A ih B, C, D, and G are defined as 
0 


(11c) 
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The closed-loop system state-space equation is then 
x = ApX + Dw (12) 

where the matrix Ac is 

Ta.-ba b s g, 

0 \ - PjM, 

The derivative matrices of A 0 . B, C, D, and G with 
respect to p arc 
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where the derivative matrices of the optimal gain matrices G s and 
Ff are given by equations (6) and (8) respectively. The derivative 
of the closed- loop state- space dynamics matrix A c with respect to 
p is 


dA c 

dp 


’dA, dB OG, d^ dG, ’ 

dp ’"dp * s ^p" # * dp 

n dA » f dK 

0 dp ’ dp M ‘ ‘ F f dp J 


(14b) 


The sensitivity of the system state and output response 3 
can be obtained from the equations 

= ~x + A-^- + (14a) 

dp dp dp dp 

dy _dr r ^x 

dp dp dp 

where the order of differentiation with respect to time t and the 
parameter p has been interchanged in (14a). 

Equations (13) and (14) can be combined into a single 
system of equations as 

B 



d£ 
dp 

or more compactly as 
x, 



dB 

dp 


y P 


= ApXp + BpW 

: ^p x p 


(15a) 

(15b) 

(16a) 
(16b) 

where the subscript p refers to sensitivity equations, and the 
definitions of the vectors and coefficient matrices follow from 
equations (15). 

Time Response - The sensitivity of the time response 
of a state-space system to known parameter variations is obtained 
by solving equations (16) as a function of time for a given input 
w(t). The theoretical solution for systems of equations of this 
type is well known 9 and is given by 

t 

(17a) 


Xp(t) = e 


V 


Xp(°) + J e 


|c Vl T> B p w(t)dx 


y p (0 = C p x p (t) 07b) 

where x p (0) is the initial condition of the system defined by 
equations (16). 

Frequency Response - The frequency response of a 
linear time- invariant state- space system can be obtained by 
Laplace transformation of the state- space equations and 
replacement of the Laplace transform variable s with the complex 
frequency s = jo) (for zero system initial conditions)^* 9 . The 
sensitivity of the frequency response can be obtained using the 
same technique on the system sensitivity equations (16). The 
result is 

y p (jto) = C p (jiol • A p ] 'BpwOto) (18) 

The frequency response h(jco) and frequency response 
sensitivity hp(j<o) of a given input/output pair calculated from 
equation (18) are complex quantities expressed in terms of real 
and imaginary components as a function of frequency ( 1 ). These 
quantities are usually more conveniently expressed in terms of 
magnitude and phase, and sensitivity of the magnitude and phase, 
than in their real and imaginary components. 


DYNAMIC RESPONSE SENSITIVITY 

Equations for the sensitivities of a given linear state-space 
dynamic system to variations in parameters which define that 
system can be obtained by partial differentiation of the state-space 
equations with respect to the desired parameter. For example, 
consider the linear, bme-invariam state-space system 
x = Ax + Bw (13a) 

y = Cx (13b) 

where x is the system state vector of dimension (n x x 1), w is the 
system input vector of dimension (n w x 1), and y is the system 
output vector of dimension (nw x 1). The matrices A, B, and C 
are appropriately dimensioned coefficient matrices. Note that 
equations (13) can describe either an open or a closed-loop 
system and that the control input vector w is taken here to be a 
general input or reference command. 


The magnitude Ihl and phase $ of a complex quantity 

h(jcu) = a + jb can be calculated from the real and imaginary 
components as 

(19a) 

<t> = tan'£ (19b) 

a 

The sensitivities of the magnitude and phase are obtained by 
differentiating equations (19) with respect to p 


aw 

dp 

d* 

dp‘ 


4L +b db 


dp 

dp 


(20a) 


(20b) 


where da/dp and db/dp are the real and imaginary components of 
dh/dp respectively. 
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Singular Values - The output vector y(s) is related to 
w(s) by a transfer function matrix of dimension n y X n u in the 
multi-input, multi-output case. When the transfer function matrix 
is some n„ x n u or n y x n y matrix H(s) of a controlled system, it 
is often desirable to compute the minimum and maximum 
singular values of H(s) as a function of the complex frequency s 

= joj, since these singular value quantities have been related to 
various control system design criterion 10 . The singular values 
and the sensitivity of the singular values to parameters as a 
function of frequency co can be obtained from H(ja>) and 

3H(jto)/dp using the definition of the singular value 
decomposition of a complex matrix. 


For the square complex matrix H of dimension nh X nh, 
the singular value decomposition of H is defined by 

H = U2V* (21) 

where (*) denotes complex conjugate transpose, U and V are 

unitary transformation matrices, and I is 

I = diagdGj o 2 ... o^]) (22) 

where di is the 1 th real scalar singular value of H. For the i lh 

singular value, equation (21) can be written as either 

Hv = u.o. (23a) 

i i i 

or 

HV = v.a (23b) 

where uj and vi are the 1 th columns of U and V respectively. 


For the case when H is a known continuous function of 
the parameter pj for which the derivative 3H/dpj is also known 
and the cq are distinct (no repeated singular values), then the 


sensitivity of the i 01 singular value 0\ to the parameter p is 11 ' 13 


dOj 


= Real 



(24) 


Eigenvalues - The eigenvalues and right eigenvectors of 
a real n x x n x matrix A are defined by the equation 9 
AE = EA (25) 

where E is the matrix whose columns are the right eigenvectors 

and A is a diagonal matrix with the i 1 * 1 diagonal element being the 
i A eigenvalue, assuming A has no repeated eigenvalues and a full 
set of linearly independent eigenvectors. The left eigenvectors 
are just the rows of the inverse modal matrix E' 1 . For the i th 
eigenvalue Xj, equation (25) can be written as either 
Ae. = X.e. (26a) 

iTa -X.il' (26b) 

where ej and 1; are the i 01 right and left eigenvectors respectively. 


The sensitivity of the I th eigenvalue Xj with respect to the 


parameter p is 14 ' 16 
3X i - - T 3A. 


dp 1 dp ‘ 


(27) 


Covariance Response • Often the response of an 
asymptotically stable linear system to random disturbances or 
inputs which can only be described in statistical terms is desired. 
In this case, the response of the system is computed using the 
covariance equation and the noise intensity matrices which model 
the random disturbance or noise 5 * 6 . Equations for the sensitivity 
of the covariance response have been developed by differentiation 
of the covariance equation with respect to the parameter variation 
of interest. 17 


The input w(t) of the linear system given by equations 
(13) is assumed to be a zero mean, Gaussian distributed, "white" 
noise with noise intensity matrix W. The steady-state covariance 
response of the (asymptotically stable) system is given by 
solution of the steady-state covariance equation 

0 = AX + XA T + BWB T (28) 

and the steady- state output covariance is 

Y = C T XC (29) 


Assuming the derivative matrices of A, B, C, and W with 
respect to a parameter pj arc known, equation (28) can be 
differentiated with respect to pi to obtain 


o= a M. + £x a t 

dp dp 


+ {^X + X 
1 dp 


3A 

dp 


+ WB T + B ^ B T - 

op dp 


(30) 


•BW^- j 
a P 

Equation (30) is linear Lyapunov equation which has a unique 
solution by virtue of strictly stable eigenvalues of the A matrix. 
The sensitivity of the output covariance is 


Q = |^xc T + 

dp dp dp 


: T + cx$£- 
dp 


(31) 


AEROSERVOELASTIC AIRCRAFT SENSITIVITY 
STUDY 


Description - An aeroservoelastic control law design 
sensitivity problem was formulated and analyzed for various 
control law and structural parameter variations. This was done to 
numerically demonstrate the application of the sensitivity 
developments of previous sections to a realistic aircraft 
stmctuiWcontrol law design problem. 


A mathematical model of the longitudinal dynamics of the 
Drone for Aerodynamic and Structural Testing, Advanced 
Research Wing-II (DAST ARW-II) Firebee aircraft. Figure 2, 
was developed for this example. 18 ' 19 The open-loop state-space 
model is of 25 th - order incorporating rigid body plunge and pitch 
modes, three elastic vibration modes, elevon and aileron control 
surfaces with actuators, and a Dryden 20 second-order vertical 
gust input model. It has elevon and aileron commanded 
deflections as inputs, pitch rate and normal acceleration at the 
c.g., wing acceleration at two locations, and actual control 
surface deflections and rates as outputs. At Mach 0.75 and 
15,000 feet altitude, the open-loop aircraft model has two real 
short period roots, one stable at -3.625/sec. and one unstable at 
1.101/sec., and a tightly damped aeroeiastic flutter mode with 
eigenvalues at -0. 162 ± j 1 1 8.3/sec. 

A control law design problem was formulated to stabilize 
the unstable short period root of the aircraft while maintaining or 
augmenting the stability of the aeroeiastic mode using reasonable 
control surface deflections and rates. Center-of- gravity pitch rate 
and acceleration, and wing acceleration from the aft wing 
accelerometer were selected as feedback measurements. These 
measurements were assumed to be noisy. A random gust 
environment of 5 ft./sec. (60 in./sec.) root mean square vertical 
gust velocity was selected for the control law design. 

The original eight mathematical model outputs were 
weighted in the LQG cost function by the matrix Q. The diagonal 
elements of Q were selected using the "Bryson" rule [reference 6, 
pg. 1691 as one over the square of the desired maximum output. 
The weighting matrix R on commanded control surface 
deflections was selected as an identity matrix. The Dryden gust 
input noise intensity was selected to give a 5 ft./sec. RMS gust 
input, and the measurement noises were selected to be 
approximately 10% of the expected output due to the gust input. 
Table 1 summarizes the numerical values for the weighting 
matrices and noise statistics. 


Numerous parameters were selected to exercise the 
sensitivity analysis techniques described in previous chapters. 
Sensitivity results for four of these parameters are presented here. 
The four parameters and their nominal values are listed in Table 
2. Parameter 1 is an element of the cost function weighting 
matrix Q, affecting the regulator portion of the optimal LQG 
control law solution. Parameters 2 is an element of the noise 
intensity matrix V, which affects the Kalman Filter portion of the 
LQG solution. Parameter 3 is a factor that simulates the effects 
of structural wing bending stiffness changes by uniformly scaling 
the two wing bending mode natural frequencies. Parameter 4 
locates the aft wing accelerometer relative to the forward 
accelerometer through a scaling of the nominal longitudinal 
separation distance between the accelerometers. Sensitivity 
matrices of the open- loop state- space model and LQG matrices to 
the four parameters were also generated. 
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Sensitivity - A numerical sensitivity analysis of the 
aeroservoelastic example problem was conducted. This 
numerical analysis consisted of the following: 1) solution of the 
optimal LQG control law problem for nominal parameter values, 
2) calculation of the sensitivity of the LQG solution to the four 
parameters, and 3) computation of the nominal controlled system 
eigenvalues, covariance response, and time and frequency 
responses, and the sensitivity of those responses, to the four 
parameters using the optimal control law sensitivity information. 

Sensitivity results presented throughout this paper are 
multiplied by the nominal value of the parameter of interest, such 
that the (semi-relative) sensitivity results for every parameter can 
be directly compared on a percent parameter change basis. This 
type of semi-relative sensitivity data has the same units as the 
response of interest in all cases. 

Closed- loop system eigenvalues and then- sensitivities to 
each of the parameters were computed using equations (26a) and 
(27). The closed'loop short period and flutter mode eigenvalues 
and their sensitivities are given in Table 3. Note in Table 3 that 
the sign order of the sensitivity of the imaginary part is 
significant. The notation ± means a positive change in the 
parameter will increase the magnitude of the imaginary part of the 
eigenvalue, whereas the notation 7+ means a positive parameter 
change will decrease the magnitude. Neither the short period nor 
flutter mode eigenvalues are affected at all by parameters 2 and 4, 
since these parameters affect only the Kalman Filter portion of the 
LQG solution. The wing bending stiffness parameter, while 
having an expected large effect on the flutter mode eigenvalues, 
also has a significant effect on the aircraft short period mode 
eigenvalues since the sensitivity results are of the same order of 
magnitude. 

Covariance responses and sensitivities of the aircraft 
model to a 5 ft./sec. RMS random vertical gust environment 
were computed using equations (28) - (31). Mean-square 
responses and sensitivities derived from the covariance results for 
aircraft pitch rate, normal c.g. acceleration, and wing acceleration 
at the forward wing accelerometer are given in Table 4. Note that 
an increase in wing bending stiffness would tend to decrease the 
aircraft pitch rate, c.g. acceleration, and the wing tip acceleration 
in the random gust environment, as would moving the aft wing 
accelerometer forward (a negative change in parameter 4). 
Parameter 1 could be used to tradeoff pitch rate response with 
c.g. and wing acceleration since the sensitivity derivatives have 
opposite signs. 

Output time responses and sensitivities were computed 
for the closed-loop aircraft subjected to the 1 - cosine discrete 
gust 

[ 60 x ( 1.0 - cosine (nt / .25)) 0.0 £ t £ .25 

w(t) = (32) 

[0.0 .25 < t < 1.0 

using equations (17). Time histories of aircraft pitch rate and 
c.g. acceleration arc shown in Figures 3 and 4. Also shown in 
these figures are the sensitivities of the responses to parameters 3 
and 4. Note here that of the two parameters, the pitch rate 
response is most sensitive to the wing bending stiffness. The 
sensitivity of the c.g. acceleration is largest with respect to 
parameter 3, indicating that wing bending stiffness is a significant 
factor in normal accelerations due to gust encounters. 


variation. These computed responses were compared with 
sensitivity derivative-based first-order predictions of the response 
computed by 

f p = fn + $ x ' ip <- v,) 

where f refers to the response of interest. Ap is the parameter 

change, 3f/3p is the appropriate sensitivity derivative, and the 
subscripts p and n refer to predicted and nominal responses 
respectively. Percentage errors in the covariance response 
predictions were calculated as 
f c -f D 

%E = -^xl00 (34) 

where the subscript c refers to the computed response. 

Validation results for the prediction of the aircraft pitch 
rate, c.g. acceleration, and aft wing accleration due to variations 
in parameter 3, the wing bending stiffness parameter, are shown 
in Figure 6. The percent errors in predicting pitch rate and c.g. 
acceleration are reasonable even for large variations in the value 
of the parameter. Further, the slope of the error curve is zero 
near the zero parameter chan ge(nominal value) point, where the 
sensitivity derivative used for the prediction was originally 
calculated. This indicates that the sensitivity derivative is exact at 
this point, verifying the derivation of the analytical sensitivity 
expressions. The percent error results for the aft wing 
acceleration prediction are larger than for the other two response 
predictions, however the error is less than about -30% for ± 10% 
variations in the parameter. In an actual application of these 
sensitivity methods, parameter variation magnitudes would 
normally be restricted by good engineering practice to be 
relatively small values, closer to the region where the sensitivity 
results are nearly exact. 

COMPUTATIONAL COSTS 


The computational burden associated with the numerical 
evaluation of the analytical LQG problem sensitivity equations 
can appear to be substantial, since solution of two matrix 
Lyapunov equations (equations (7) and (9)) is required to obtain 
the sensitivities of the linear quadratic regulator and Kalman Filter 
gain matrices to a single parameter. For this reason, a 
comparison of the analytical sensitivity evaluation versus onc- 
and two-step finite difference calculations for the equivalent 
sensitivity information was made. The measure of comparison 
for the three calculations was central processing unit time (CPU 
seconds) on a Digital Equipment Corporation MicroVAX II 
computer, where the LQG sensitivity equations were 
programmed as user functions to a commercially available linear 
systems analysis computer code. 


The one-step regulator gain matrix G finite difference 
sensitivity was calculated as 


8G _ °1 G n 
8p ' .025 


(35) 


where Gj is the perturbed LQR gain for the wing bending 
stiffness parameter (parameter 3) perturbed positively by 2.5% 
(L025 times the nominal value). The KF gain matrix F finite 
difference sensitivity was calculated similiarly. 


The complex frequency response and sensitivities of the 
elevon open -loop transfer function with the aileron loop closed 
were calculated using equation (18). The complex (real and 
imaginary) results were converted to magnitude and sensitivity of 
the magnitude using equations (20). The magnitude result is 
shown in Figure 5, as arc the sensitivities of the magnitude to 
parameters 3 and 4. The magnitude is most sensitive to the wing 
bending stiffness at about .6 rad/sec., although the peak 
sensitivity for parameter 4 coincides with the peak of the 
magnitude at 1. 1 rad/sec. 

Sensitivity Validation - The sensitivity analysis 
results were evaluated for accuracy by comparing predicted 
covariance responses with covariance responses computed for 
variations in the nominal values of the parameters. The four 
parameters were varied ± 25% from their nominal values in 5% 
increments, and the new optimal control law and controlled 
aircraft covariance response were computed for each parameter 


The two-step regulator gain matrix G finite difference 
sensitivity was calculated as 


80 _ G 2' G i 
8p " .05 


(36) 


and the wing bending stiffness parameter was perturbed by ± 
2.5%. G 2 refers the LQR gain matrix obtained for parameter 3 
perturbed to 1.025 times nominal, and Gi refers to the gain 
matrix for parameter 3 perturbed to .975 times nominal. The 
Kalman Filter gain matrix F sensitivity was again calculated in the 
same manner. 


The results of the CPU time comparisons are shown in 
Table 5, as is the CPU time required for solution of the LQG 
problem without sensitivity calculations. These results show that 
the analytical sensitivity expressions require substantially less 
CPU time than either the one- or two-step numerical finite 
difference approaches for a single parameter sensitivitv analysis. 
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and require only a 33% increase in CPU time over the nominal 
LQG problem solution. Furthermore, the computational 
advantage of the analytical approach is likely to increase when 
sensitivity calculations for more than one parameter are involved, 
since additional computational efficiency can be achieved by 
storage of the decomposed coefficient matrices of the sensitivity 
equations (7) and (9). Similar computational efficiencies are not 
possible with the finite difference approaches, since they require 
a solution of the LQG problem for each perturbation of the 
parameter of interest 

CONCLUSION 

This paper has highlighted a method for computing the 
sensitivity of optimal LQG control laws to various parameters 
using analytical sensitivity expressions. The LQG sensitivity 
results are used to predict changes in closed-loop aircraft 
responses due to changes in the nominal values of the parameters 
of interest These sensitivity results are shown to be useable for 
hierarchal integrated structure/control law design problems 
through a large aeroservoelastic aircraft example. Sensitivities of 
covariance, time, and frequency responses of the aircraft to 
various parameters were computed. The sensitivity results were 
validated against computed response changes due to changes in 
the nominal values of various parameters and found to be 
accurate for ± 15% changes in the parameter values. It was also 
found that it is cheaper to evaluate the analytical LQG sensitivity 
expressions than to calculate the equivalent sensitivity 
information by finite difference means. 
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TABLES 


Table 1 LQG Control Law Problem Data 


Svmbol 

Value 

Description 


"1TDT ' 

Pitch Rate Weight 

0(2.2) 

1.00 

C.G. Acceleration Weight 

0(3.3) 

0.01 

Fwd. Wing Acceleration Weight 

0(4,4) 

0.01 

Aft Wing Acceleration Weight 

0(5,5) 

2.04 x 10-' 

Elevon Deflection Weight 

0(6.6) 

4.40 x 10‘ 3 

Aileron Deflection Weight 

0(7.7) 

1.56 x 10- 4 

Elevon Rate Weight 

0(8.8) 

1.83 x 10-6 

Aileron Rate Weight 

R(l.l) 

1.00 

Commanded Elevon Weight 

R(2,2) 

LOO 

Commanded Aileron Weight 

V(l.l) 

2.00 x 10" 3 

Pitch Rate Noise Intensity 

V(2,2) 

1.50 x lO' 3 

C.G. Accl. Noise Intensity 

V<3,3) 

6.00 x 10- 3 

Aft Wing Accl. Noise Intensity. 

U(l,l) 

1.00 x lO 6 

Input Noise Intensity (Elevon Loop) 

U(2,2) 

1.00 x !0- 6 

Input Noise Intensity (Aileron Loop) 

W 

3.60 x 10 3 

Dryden Gust Model Noise Intensity 

Table 2 

Integrated Aircraft Design Problem Parameters. 

Number 

Nominal Value 

Description 

1 

0.01 

Q( 1,1) Pitch Rate Weight 

2 

2.00 x iO* 3 

V(l,l) Pitch Rate Noise Intensity 

3 

LOO 

Wing Bending Stiffness Factor 

4 

7.58 

Aft Wing Accl. Location 

Table 3 

Semi-Relative Closed-Loop Short Period and Flutter 


Mode Eigenvalue Sensitivities to Parameters. 


Short Period Mode Flutter Mode 


(1/sec.) 

(1/sec.) 

Param. 

(-5. 1 36 + j2.742) (-5.046 ± j 1 . 178 x 10 2 ) 

1 -4.74 x 10-1 ./+ J7.28 x 10-' 4.69 x lO" 5 -/+ jl.50 x 10*5 

2 

0.00 

0.00 

3 

4.72 x 10-' ±j4.04 3.70 ± j2.83 x 10' 

4 

0.00 

0.00 


6 


Tabic 4 Semi-Relative Gosed-Loop Mean- Square Response 


Sensitivities to Parameters (5 ft/scc. RMS Gust). 


Param. 

Pitch Rate 

(deg./sec.) 
(5.15 x 10- 2 ) 

C.G. Accl. 
(g's) 

(2.65 x 10-2) 

Wing Accl. 
(g's) 

(2.35 x 10>) 

1 

-6.18 x 10*3 

2.24 x 10" 4 

6.24 x 10- 5 

2 

5.91 x 10- 3 

9.91 x Ifr 4 

5.21 x 10-3 

3 

-9.84 x 10-2 

-3.38 x 10-3 

-5.32 x 10' 

4 

1.35 x 10- 3 

1.16 x 10-3 

2.94 x 10-' 


Table 5 


CPU Time Comparisons of Analytical Sensitivity 

TuneTSecT 

LQ<j Solution Only 


Analytical Sensitivity Expressions 
One-Step Finite Difference 
Two-Step Finite Difference 


m.6r 

133.55 

196.44 

287.39 


Pitch Rate 

q 




























t I 1 1 

0.0 0.25 0.5 0.75 1 .0 


Time t (sec) 



0.01 0.1 1.0 10.0 100.0 1000.0 


Frequency o> (rad/sec) 


Figure 4. C. G. acceleration response to 1 - cosine discrete 
gust and sensitivity to parameters 3 and 4. 


Figure 5. Open elevon loop transfer function magnitude and 
sensitivity to parameters 3 and 4 (aileron loop 
closed). 
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% Change Parameter 3 



% Change Parameter 3 


% Error in 
Predicted Aft 
Wing Acceleration 
Response 



Figure 6. Percent prediction errors in covariance responses for 
± 25 % changes in parameter 3 nominal value. 
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